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While the particle hypothesis for dark matter may be very soon investigated at the LHC, and as 
the PAMELA and GLAST satellites are currently taking new data on charged and gamma cosmic 
rays, the need of controlling the theoretical uncertainties affecting the possible indirect signatures of 
dark matter annihilation is of paramount importance. The uncertainties which originate from the 
dark matter distribution are difficult to estimate because current astrophysical observations provide 
rather weak dynamical constraints, and because, according to cosmological N-body simulations, 
dark matter is neither smoothly nor spherically distributed in galactic halos. Some previous studies 
made use of N-body simulations to compute the 7-ray flux from dark matter annihilation, but such 
a work has never been performed for the antimatter (positron and antiproton) primary fluxes, for 
which transport processes complicate the calculations. We take advantage of the galaxy-like 3D 
dark matter map extracted from the HORIZON Project results to calculate the positron and an- 
tiproton fluxes from dark matter annihilation, in a model-independent approach as well as for dark 
matter particle benchmarks relevant at the LHC scale (from supersymmetric and extra-dimensional 
theories) . We flnd that the flux uncertainties arise mainly from fluctuations of the local dark matter 
density, and are of ~ 1 order of magnitude. We compare our results to analytic descriptions of the 
dark matter halo, showing how the latter can well reproduce the former. The overal antimatter 
predictions associated with our benchmark models are shown to lie far below the existing measure- 
ments, and in particular that of the positron fraction recently reported by PAMELA, and far below 
the background predictions as well. Finally, we stress the limits of the use of an N-body framework 
in this context. 
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I. INTRODUCTION 

The idea that dark matter is made of exotic weakly 
interacting massive particles (WIMP) is very appealing in 
the sense that the related existing frameworks in particle 
physics beyond the standard model (BSM) could solve (i) 
problems inherent to elementary particle theory, e.g. the 
force unification scheme and/or the stabilization of the 
theory against quantum corrections, which would thus 
confer a more fundamental meaning to this theory; and in 
the sametime (ii) the dark matter issue as characterized 
in astrophysics and cosmology (see e.g. [l[ for a recent 
review). Indeed, the cold dark matter (CDM) paradigm 
seems to feature a powerful and self-consistent theory of 
structure formation, which at present is able to reproduce 
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and explain most of large and mid scale observations (see 
e.g. and references therein). 

While the Large Hadron Collider at CERN (LHC) is 
about to start hunting BSM particle physics and may 
soon provide new insights about the dark matter parti- 
cle hypothesis, high energy astrophysics experiments are 
also of strong interest to try to determine the actual na- 
ture of dark matter. Indeed, if dark matter is made of 
self-annihilating particles, a property which is found in 
many BSM models and which provides a natural mecha- 
nism to explain the dark matter cosmological abundance 
as observed today, there should be traces of annihilations 
imprinting the cosmic ray spectra (indirect detection of 
dark matter, see e.g. II, i, i, Hi)- Therefore, sateUite 
experiments such as PAMELA and GLAST, which are 
dedicated to large field of view observations of charged 
cosmic rays and 7-rays, respectively, in the MeV-TeV en- 
ergy range, should be able to yield additional constraints 
or smoking guns very soon in this research field 0] ■ The 
most promising astrophysical messengers that could trace 
the annihilation processes are indeed gamma-rays and 
antimatter cosmic rays, which have been first investi- 
gated in 0] and in Q respectively. Confinement and 
diffusion on magnetic turbulences limit the origin of the 
latter to our Galaxy, whereas the former, which only ex- 
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periences the common r ^ flux dilution, may be observed 
even when emitted from extra-galactic regions. 

Nevertheless, predictions of these exotic signals are af- 
fected by many uncertainties coming from (i) the under- 
lying WIMP model (ii) the distribution of dark matter 
in the relevant sources (iii) the propagation of charged 
cosmic rays in the Galaxy in the case of searches in the 
antimatter cosmic ray spectra. The first point is encoded 
in the WIMP mass, the annihilation cross-section and 
the annihilation final states which fully define the in- 
jected cosmic ray spectra: this has been widely inves- 
tigated for years. The second point is currently sur- 
veyed by the state-of-the-art N-body experiments, but 
the use of N-body dark matter maps in the context of 
indirect detection has only been performed for 7-ray pre- 
dictions [13, El, [H, Besides, there has been lots of 
efforts to estimate the effect of varying the dark matter 
distribution by means of analytical calculations, allow- 
ing for instance to study the effects of sub-halos on the 
7-ray [3, fisl [igI. ITtI . [ist as well as on the antimatter pri- 
mary fluxes [i^, [2^ . Such analytical studies mostly rely 
on the statistical information supplied by N-body simu- 
lations. However, though they provide very nice theoret- 
ical bases to unveil the salient effects of any change in 
the relevant parameter space, and permit to go beyond 
the current numerical resolution limits, they are usually 
bound to simplifying hypotheses, such as the sphericity 
of sources. Especially, no study of the antimatter signa- 
tures has been directly performed in the frame of N-body 
environments up to now. Finally, point (iii), referring to 
cosmic ray propagation, has also widely been investigated 
for antiprotons [2l|, H^l and positrons [13, H^l , and while 
the transport processes depend on sets of still degener- 
ate parameters within different propagation models, the 
origins of uncertainties are now rather well understood. 

In this paper, we make use of a realistic (in the sense 
of ACDM cosmology) 3D dark matter distribution of a 
galactic-sized halo extracted from the HORIZON Project 
simulation results to study the antimatter signals. Our 
main purposes are to characterize and quantify the uncer- 
tainties coming from dark matter inhomogcneitics (not 
necessarily clumps) and departure from spherical sym- 
metry in a virialized Milky- Way-like object; and finally 
make predictions for different dark matter particle can- 
didates. Despite the rather low resolution of our numeri- 
cal data compared with the highest-level artillery on the 
market, as portrayed, e.g., by the Via Lactea simula- 
tions [HI, UM I study for the first time the production 
of antimatter cosmic rays from dark matter annihilation 
in an N-body framework. The resolution issue is however 
less important for charged cosmic rays than for gamma- 
rays, because the signal is strongly diluted by diffusion 
effects, which therefore tends to smooth the yields from 
high density fluctuations, provided the latter are not too 
close to the observer. Nevertheless, as very local density 
fluctuations are expected to affect the antimatter signals, 
we will quantify and correct the consequence of loss of 
resolution by analytically extrapolating our results. For 



the WIMP candidates, we will use a model-independent 
approach as well as typical models of BSM particle the- 
ories, some of them being observable at the LHC. 

This article is parted as follows: we recall the salient 
features of the HORIZON N-body experiment in Sect. [Til 
then, in Sect. IIIIl we briefly review the positron and 
antiproton propagation before sketching the method we 
adopt to connect the propagation to the N-body source 
terms; we describe the WIMP models in Sect. IIVI we fi- 
nally present and discuss our results and predictions of 
the primary positron and antiproton fluxes in Sect. |V| 
(where the expert reader is invited to go directly) before 
concluding. 



II. THE HORIZON FRAMEWORK 

The framework of this study is a simulation of the Hori- 
zon collaboration. The run used the Adaptive Mesh Re- 
finement code RAMSES [2^ with an effective number of 
particles of Np ~ 1024^ in a box of size L = 20/i~^ Mpc 
where the initial conditions were fixed by the WMAP3 
results (17„j = 0.24, f^A = 0.76, 0^ = 0.042, n ^ 0.958, 
Ho = 73, as = 0.77). 

A Milky- Way sized halo at 2 = was selected and re- 
fined using the so-called "zoom" technique. We obtained 
a maximum linear resolution of about 200 pc and the 
particle mass is Mp = 7.46 10^ in solar mass (Mq) units 
(we refer the reader to ref. [l^l for more details on the 
characteristics and the analysis of the simulation.). 

A non biased calculation of the dark matter density 
around a simulation point i is given by the algorithm of 
Casertano and Hut [26| . namely: 



(1) 



where V{rj) = 47r/3r^ is the volume of the smallest 
sphere around the particle i that includes j neighbors. 

The virial radius (defined as p{ryir) = 200pc) of our 
halo is equal to 253 kpc, corresponding to an enclosed 
mass of 6.05 x 10"Mq or 8.1 x 10^ particles. Within this 
radius, the dark matter halo density can be fitted by the 
usual spherical parameterization: 



P{r) = Po 



1 + [r/af 



(2) 



where po = Po is the local density in the solar neigh- 
borhood and To = i?0 = 8 kpc is the distance from the 
Sun to the Galactic center (GC). Density maps, after 
projections in the {x,y) and {z,x) planes, are shown on 
Fig. [1] where the Earth has been located at three differ- 
ent positions: at x = 8, y = 8 and 2 = 8 kpc, respec- 
tively (see Sect. IIII Bl for further details). When consid- 
ering the radial distribution of our data, the best fit gives 
(a,/3,7) = (0.39,3.72,0.254), and a = 13.16 kpc, but a 
NFW-like profile (1,3,1) - (a = 10 kpc) and a cored 
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profile (0.5,3.3,0) — (a = 4.5 kpc) are also in good sta- 
tistical agreement due to the resolution limits, especially 
in the central region. The averaged density at a radius 
of 8 kpc is pq = 0.25 ± 0.19 GeV/cm^. 108 sub-halos 
have been identified with masses between 2 x IO^A/q 
and 2.4 x lO^^M© with a mass distribution following 
N{Msub) oc MJ^^ above 5 x IO^Mq. Inside those sub- 
structures, the dark matter density scales like a universal 
power law Pdii") r~^'^ for the outer part. Regarding 
the inner part, the logarithmic slope is highly specula- 
tive due to the resolution limit. Furthermore, all clumps 
are found to be rather far from the Galactic center and 
will actually not modify our coming results. Anyway, 
though we will further mimic higher resolutions thanks 
to analytical extrapolations, this N-body environment is 



still relevant for the punch line of this study, which is 
to quantify the importance of solar neighborhood dark 
matter density fluctuations. 



III. ANTIMATTER COSMIC RAYS 

Antimatter cosmic rays are very interesting messengers 
to look for dark matter annihilation traces, because they 
arc seldom produced in standard astrophysical processes. 
The astrophysical secondary background at the Earth is 
generically created out of the spallation of cosmic rays 
oflF the interstellar medium (ISM), and is expected to be 
highly power-law suppressed at energies above ~ 10 GeV. 



For charged Galactic cosmic rays in the GeV-TeV regime, the relevant (steady state) propagation equation which 
characterizes the transport, for any species, of the cosmic ray number density per unit of energy Afcr = dnc-[{E) / dE 
reads: 



K{E)VN„ - FconvAAer 



-I- 



d_ 
'dE 



b{E)K. 



d 

oE 



T{E)J\f„ + Q = ; 



(3) 



where K{E) encodes the spatial diffusion due to scattering off magnetic turbulences, T4onv is the velocity of the 
convective wind that drifts cosmic rays out of the Galactic disk, b{E) stands for energy losses, Kee characterizes 
diffusion in momentum space, r(i?) features the spallation reactions and Q is the (stationary) cosmic ray source of 
interest. I 



Though the above equation stands for any cosmic ray 
species, some of the processes can be safely neglected in 
the GcV regime, depending on the species. For positrons, 
the main processes are spatial diffusion and energy losses, 
so that we will neglect convection and reacceleration. 
Such an approximation is very often used in the litera- 
ture, and its relevance is explained in great details in [2^ . 
For antiprotons, the dominant processes beside diffusion 
are convection and spallation, which are more efficient at 
low energies, typically < 10 GeV, whereas reacceleration 
and energy losses do not play a major role [2l|, [2I]. We 
will therefore neglect those latter effects for antiprotons 
in the following. 



otic flux at the Earth is merely given by: 

#cr(-E) ^ Vcr_ „ 

dE ~ Att ^ 

X I d^xsG{E; xq ^ xs) X Q{xs) ,(4) 

J slab 

where Vcr is the cosmic ray velocity, S = {{<J annv) / 'i) x 
(Pq/ttt,^)^ allows a convenient normalization to the lo- 
cal WIMP-annihilation rate, Q is the spatial source 
term that we will further discuss in Sect. IIII Bi and 
Q = J dEsQdN/dEs is the Green function convoluted 
with the injected cosmic ray spectrum dN/dEs- 

In the following, wc present the main features of our 
propagation model for both antiprotons and positrons. 
We will end this section by discussing the source term Q 
as set in our N-body framework. 



The diffusion equation can be solved cither fully nu- 
merically (see e.g. [Hf) or semi-analytically in typical 
energy regimes and/or geometries of the diffusion zone 
(see e.g. [301). Wc will use the latter approach in this pa- 
per. Generic solutions for point-like sources oc 6'^{Es, xs) 
(4-D Dirac function) can be expressed in terms of Green 
functions G{E, Xq ^ Es, xs) that characterize the prob- 
ability for cosmic rays of energy Es injected at position 
Xs) to be detected at the Earth with energy E. Hence, 
for a WIMP particle of mass and of thermally aver- 
aged annihilation cross-section ((Jannw), the primary ex- 



A. Propagation model 

The propagation model that we adopt in this paper 
is characterized by a slab geometry with the origin lo- 
cated at the Galactic center, a radial extent i?siab and 
a vertical half-height L. This defines the cosmic ray 
confinement zone (see Fig. [T]), thanks to which we im- 
pose Dirichlet boundary conditions to Eq.([3]), that is 
A/'(-Rsiab, z) = 7V(r, L) = 0. The Earth is located on a 
circle in the Galactic disc with cylindrical coordinates 
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X axis [kpc] - N-body frame y axis [i<pc] - N-body frame z axis [kpc] - N-body frame 



Figure 1: Projected dark matter density in the cosmic ray diffusion piane, for different choices of the Earth iocation (Cartesian 
coordinates). Left: Earth at x = 8 kpc (p© = 0.23 GeV.cm"'^), diffusion piane in Oxy; middie: Earth at y = 8 kpc 
(p0 = 0.28 GeV.cm"''), diffusion plane in Oxy; right: Earth at z = 8 kpc (pq = 0.13 GeV.cm^'^), diffusion zone in Oxz. The 
diffusion slab contour is explicitly drawn (visually ellipsoidal instead of spherical due to different horizontal and vertical scales) . 
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Table I: Propagation parameters compatible with B/C data 
giving the median antiproton primary (from dark matter an- 
nihilation) fluxes ; from [2H . 



(r, z) = (8,0) kpc (the polar angle is left free for the 
moment). As regards the transport processes, we take 
a spatial independent diffusion coefficient of the form 
K{E) = PKoH^ (where /3 and 7^ = pc/Ze are the parti- 
cle velocity and rigidity respectively) and a constant wind 
14onv directed outwards the Galactic plane along the ver- 
tical z axis. The reader is referred to [3l| for a detailed 
presentation of this framework. The propagation param- 
eters can be constrained with cosmic ray measurements 
of secondary to primary ratios, and are found to be de- 
generate [3l|, [22, [331 • Throughout this paper, we will use 
the median set of propagation parameters defined in [2l] | , 
which has been constrained with the current B/C data 
and which is recalled in Tab. [D Once the fluxes will be 
computed, they will be modulated to account for the so- 
lar modulation. We will use the force field approximation 
method, with an electric potential (f> = 600 MV (33 |. 



1. Antiproton and positron propagation 



The exotic antiproton production an d p ropagation 
have widely been studied in the literature 3^j_ 3^^ 3^, 



1381, and this also holds for positrons 0, [2a[3a j^El ! 
In this paper, we will use the Green functions as de- 
scribed in [2(| for both the propagations of positrons and 
antiprotons. Both species experience the same diffusion 



processes, but have quite different propagation behav- 
iors: positron propagation is mainly affected by Inverse 
Compton (IC) energy losses, while antiprotons, though 
their energy losses are negligible, are swept out of the 
galactic disc (destroyed) by the convective wind (spalla- 
tions off the interstellar gas) at low energy. Therefore, 
their respective horizon is completely reversed one from 
the other in terms of the energy dependence. 

The typical propagation length for positrons is given 
by [H: 



1/2 



(5) 



where te — lO^^s is the characteristic timescale of the 
energy loss in the IC regime, e = (E/lGeV) stands for 
the energy at which the positron is detected at the Earth, 
while es > e is the energy at which the positron is in- 
jected at the source. Given £5, this is clearly a decreasing 
function of the energy, which is hardly greater than a few 
kpc. 

For antiproton, convection is the relevant quantity that 
bounds the propagation length, and the horizon is set 
by El: 



A 



D 



K{E) 

^onv 



(6) 



which is. contrarily to the positron case, an increasing 
function of the energy, and which can quickly get values 
much greater than the typical size of the diffusion zone. 

The characteristic lengths for both antiprotons and 
positrons are reported in Fig. [21 as functions of the de- 
tected energy at the Earth. For positrons, we have con- 
sidered different injected energies Es of 250, 500 and 1000 
GeV, and A/? is plotted as a function of the fraction of 
Eg still carried by the particle at the Earth. For antipro- 
tons, we have plotted A^) as a function of the detected 
energy divided by 1 TeV. 



5 



Characteristic diffusion length (for L = 4 l<pc) 


anti-proton [E 


= 1 TeV] 


positron [E^ = 


1 TeV] 


positron [Eg = 


0.5 TeV] 


positron [E = 


0.25 TeV] 




1 1 ' 



10^ 10^ 10^ 11 



[E /E„fore*]or[E/E forp] 

D S ' max ^ 

Figure 2: Propagation scales for positrons and antiprotons 
as functions of energy. For the former, the energy reported 
on the X axis is Ed/Es, that is the detected energy divided 
by the injected energy (positrons loose energy), and for the 
latter, this is merely E / Ema.yi- Scales are normalized to L = 4 
kpc, the vertical half-height of the diffusion zone. 



These quantities are very useful to figure out the vol- 
ume from the Earth we are sensitive to when computing 
the exotic antimatter fluxes. When the positron (an- 
tiproton, respectively) energy is high (low), contributions 
to the flux are mostly local, and large fluctuations are 
therefore expected if the Earth position is varied inside 
an N-body environment because the dark matter den- 
sity varies. This is the contrary at low (high) energies, 
for which contributions are integrated inside much larger 
volumes, and tend to smooth the density gradients. 

This will be discussed in more details in Sect. IV A 21 
and IV A3| but the reader can already look at Fig. IH 
where it is shown that different local distributions of dark 
matter strongly affect the positron and antiproton sig- 
nals, at high and low energy respectively. 



B. Connection to the N-body source term 

N-body simulation results consist in a set of points 
with coordinates and masses. Though all points corre- 
spond to particles of the same mass, the algorithm of 
density computation dcflncd in Sect, [ill allows to asso- 
ciate a dark matter density to each point, which actually 
reflects the particle number density in the point neigh- 
borhood. Therefore, each point will be considered as an 
individual source term of Eq. ([3|). In order to compute 
the antimatter flux at the Earth, we have to fix the lo- 



cations of the Earth and of the diffusion zone with re- 
spect to the N-body system of coordinates. Due to the 
cosmological origin our simulation, we have no freedom 
in rescaling the distances, so that the Earth is set at 
i?0 = i?0 = 8 kpc, while the Euler angles are random 
parameters. Once the Earth position is set, we also have 
to randomly fix the slab diffusion zone, where obviously, 
the Earth itself is located. 

As already detailed in [23|, constraints on the local 
dark matter density are very important for indirect de- 
tection of dark matter with antimatter cosmic rays, be- 
cause local contributions dominate the signal in energy 
ranges depending on the species (high/low energy for 
positron/antiproton). Therefore, it is meaningful to add 
a prescription on top of the random procedure described 
above, so as to account for the observational constraints 
on the local dark matter density. According to [lOl (and 
references therein), it is reasonable to bracket pq within 
[0.2,0.4] GeV.cm"^. Notice that such a constraint is by 
itself a very strong limit on the possible variations of the 
expected fluxes. This wfll be discussed in Sect. IVl 

For visual purpose, we show in Fig. [1] projections of 
the dark matter density p in the diffusion slab according 
to three different positions for the Earth, at 8 kpc along 
each axis x, y and z, respectively (the associated local 
densities are = 0.23, 0.28, 0.13 GeV/cm^). These 
can be seen as different maps of the source term of the 
diffusion equation ^ (the actual source term is in fact 
more contrasted because scaling like p^). Note that those 
three different configurations arc typified by very differ- 
ent local environments, with sizable density fluctuations, 
and one can expect some variations on the predicted an- 
timatter cosmic ray fluxes (to go straight to the results, 
see Fig. m. 



IV. WIMP MODELS 

Despite its numerous successes, the high energy physics 
standard model of elementary particles needs some ex- 
tensions. Indeed, some questions like the hierarchy prob- 
lem, the force uniflcation or the neutrino masses can not 
be answered without advocating new physics. Further- 
more in a large variety of standard model extensions, new 
particles should be present and some of them could be 
WIMPs and candidates to dark matter. Some of those 
scenarios have deep theoretical motivations. Supersym- 
mctric models like the Minimal Supersymmetric Stan- 
dard Model (MSSM) are inspired from supergravity or 
string theory and predict new particles below the TeV(s) 
scale. Those theoretical frameworks also require the ex- 
istence of extra spatial dimensions, and have inspired the 
building of phenomenological models in which new parti- 
cles related to those new dimensions can arise at energies 
as low as the TeV scale. Some more effective models 
propose to extend the standard model in more minimal 
ways, in order to address some specific questions and 
to consider new phenomcnologies without hypotheses on 



6 



the tricky and unsolved question of the hnk with some 
possible fundamental aspects or inspiring sector. After a 
model independent setting, we will describe some typical 
benchmark situations of the MSSM, Kaluza-Klein dark 
matter, Little Higgs Model and Inert Doublet Model, of 
which we will further survey the signatures in the anti- 
matter cosmic ray spectra. 

Note that all cosmic ray spectra for the WIMP mod- 
els discussed in this paper have been generated with 
the PYTHIA Monte Carlo and some of the asso- 
ciated relic densities have been computed with the mi- 
crOMEGAs code [3]. 

A. Model-independent setting 

To allow discussion in a very general context, we define 
here two WIMP fiducial models associated with positron 
and antiproton production, respectively. For both mod- 
els, we will take a WIMP mass of 200 GeV and a total 
annihilation cross section of {(Jannv) = 3 x 10~^^ cm^/s. 

The first one, associated with positrons, is either 
a bosonic or a Dirac fermionic WIMP and is sup- 
posed to provide monochromatic electrons and positrons 
by fully annihilating like XX ~^ e+e^. This might, 
though very roughly, mimic Kaluza-Klein dark matter 
(see Sect. IIVB 21) . Such a model greatly simplifies the 
calculation for the positron flux since the source term is 
(X d{E — m-)^), and is very convenient to discuss the main 
features of our results. 

The second one, associated with antiprotons, is a 
Majorana fermionic WIMP that fully annihilates in bb. 
Such a model is generically found in SUSY theories (see 
Sect. IIVB II) . 

B. Models observable at the LHC 

We will show the resulting cosmic ray fluxes from 
WIMP annihilation in some popular BSM frameworks 
relevant at the LHC scale that we describe in the next 
paragraphs. 

1. The Minimal Supersymmetric Standard Model 

The most popular new symmetry of high energy 
physics is the supersymmetry (SUSY) [l^. The simplest 
extension in this way is known as the Minimal Supersym- 
metric Standard Model (MSSM) The gauge group is 
still SU (3) X SU{2) xU{l) and the supersymmetric trans- 
formations associate a scalar (fermion) to each fermion 
(boson) of the standard model. This model unifies the 
coupling constants of the three forces at high energy 
and nicely achieves the electroweak symmetry breaking 
through radiative corrections. It requires two Higgs dou- 
blets and we note fi the mass parameter of those (su- 
pcr)ficlds and tan /3 the ratio of their vacuum expected 



value (vcv). This model has numerous (106) parameters 
but is usually simplified assuming at high energy unified 
values for scalar (mo) and fermions (TO1/2) masses and for 
trilinear couplings (Aq). Thanks to the requirement of ra- 
diative electroweak symmetry breaking, the last param- 
eter is the sign of fi. The set {mp, tf^n/?! 5ign(/^)} 
defines a so called Constrained Minimal Supersymmetric 
Standard Model [13, IS (CMSSM or mSUGRA) point. 
In the main part of this parameter space, the dark mat- 
ter particle is the Hghtest neutralino, xi = X = DM. 
The four neutralinos Xi,2,3,4) are the mass eigenstates 
coming from the mixing of neutral Higgs and gauge bo- 
son superpartners. There were many studies dedicated to 
predictions for antimatter fluxes with SUSY dark matter, 
see e.g. [Ulilliilll. 

Inspired from the SnowMass Point and Slope [5l| defi- 
nition, we will study in the following some typical bench- 
mark points. We have chosen three points as close as 
possible to SnowMass points but respecting the WMAP 
constraint. For the sake of completeness, we also re- 
laxed the GUT scale universality assumptions to obtain 
two additional points providing light neutralinos inspired 
by jUjIll] and potentially observable at the LHC (sec the 
main parameters in Tab. [11)) . 



2. Extra dimensions 

Initially inspired by the hierarchy problem, large ex- 
tra dimensions with testable (< TeV) phenomenologi- 
cal consequences have been studied during the last ten 
years. For dark matter phenomenology, two kinds of 
models have been proposed. Either universal large com- 
pact extra dimensions with flat geometry [s^l where all 
the Standard Model fields may propaga te, or extra di- 
mensions with warped geometry [551 . |56| where all fields 
but the Higgs are in the bulk. 

New particles arise as the Kaluza-Klein (KK) modes 
associated to each standard field. 

In universal extra dimensions (UED), the KK parity 
(coming from translational invariancc) implies that the 
lightest KK particle (LKP) is stable. Consequently, the 
first mode of neutral particles (mainly photon and Z but 
also neutrino) can be interesting WIMP candidates. The 
tree level masses are degenerate and = n/R, where R is 
the size of the compact dimension and n is the mode num- 
ber. The radiative corrections add extra terms breaking 
the degeneracies and in the gauge boson sector, the elec- 
troweak symmetry breaking induces mixing between the 
gauge eigenstates. These contributions lead to different 
masses for the Z and photons modes. The dark matter 
particle is most likely a gauge boson. Here, we will con- 
sider the LKP dark matter model of [s^l, which is the 
first KK excitation of the B^^^ boson. Previous analyses 
involving this model in the frame of antimatter cosmic 
rays can be found in e.g. [s^, Isol. Issj. 

In warped geometry, one has to impose a Z3 symmetry 
to preserve proton stability. Consequently, the lightest 
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Zj, charged particle (LZP) is also stable and has to be 
the KK RH neutrino to be a WIMP candidate [li,!!^]. 
The dark matter particle is thus a Dirac fermion. The 
4-D masses arise from the localization of wave function of 
massless modes along the S*'' dimension. This mechanism 
is popular to explain the structure of Yukawa couplings. 
For previous studies of indirect detection involving LZP 
DM candidates, see e.g. [la. .37. .39. 61]. 

We have selected one UED model and another one in 
warped geometry (see Tab. HI)). 

3. The Inert Doublet Model 

The scalar sector is a key point of numerous high en- 
ergy physics theories, especially for the standard model 
and some of its extensions. Up to now, no scalar particle 
has been discovered and this sector is still unknown. The 
Inert Doublet Model [H^l is a more ad-hoc framework, 
but quite minimal. It is a two Higgs doublet extension 
of the standard model with a Z2 symmetry under which 
one of the doublet and the fields of the Standard Model 
are even. The Zi symmetry is not broken preserving 
flavor changing neutral currents and forbidding a vac- 
uum expectation value for the odd doublet. This scalar 
sector provides three additional odd particles, a neutral 
{Hq = H = DM), a pseudo (Aq) and a charged scalar 
In our notation, a point of this model is then de- 
fined by the mass of the extra scalar particle masses : 
Mho j , Mh+ , the mass of the standard model Higgs 
rrih and two scalar potential parameters related to the 
odd doublet: a mass dimensional parameter fi2 and a 
coupling A2 (see (63j for more details concerning the pa- 
rameterization) . The neutral extra scalar particle of the 
model, Hq, is a typical candidate of WIMP dark matter 
with an interesting phenomenology. Inspired from an ex- 
tended study of the dark matter phenomenology of the 



Inert Doublet Model [6^ we have taken two benchmark 
points respecting the WMAP constraint (see Tab. HI]). 

4- Little Higgs Model 

Another possibility to solve the hierarchy problem is 
to build a model such that the Higgs particle appears 
as a Nambu-Goldstone boson. In the Little Higgs model 
with T-parity (see [13, [sl, [1^ and references therein), 
a global symmetry SU(5) with locally gauged subgroup 
[SU{2) X U{1)]'^ is spontaneously broken down to S0{5). 
The local subgroup is also broken down to the diag- 
onal subgroup SU{2) x U{1), identified with the SM 
gauge group. To prevent the contributions from the new 
heavy gauge bosons associated with the enlarged sym- 
metry to spoil the electroweak precision measurements, 
a Z2 symmetry, the T-parity, is introduced (67j . There- 
fore, the lightest T-odd particle, which happens to be the 
heavy photon Ah , is a potential candidate for dark mat- 
ter [6^. Other T-odd heavy bosons include the heavy 
gauge bosons Wh and Zh and a scalar triplet Higgs bo- 
son $. The parameters of a point in this model will 
be given by the value of the vcv associated with the 
symmetry breaking /, the masses of the heavy T-odd 
bosons, itiahj ttt-Wh^ '^Zh ^'^d m$, and the mass of the 
Higgs particle rrih ■ The annihilation of the heavy photon 
leads to the production of W , Z or h. Points that satisfy 
the WMAP constraint lie on two branches in the plane 
/ — m/i. For points on the upper branch (U-branch), 
rrih > 2m Ah , while for points on the lower branch (L- 
branch), rrih < 2m a„ (69[. Therefore, we have selected 
one benchmark point on each branch (see Tab. HI)) . 

To summarize our benchmarks, we have chosen five su- 
persymmetric points, two models with extra dimension, 
two Inert Doublet points and two Little Higgs models. 
Parameters and details can be found in Tab. [IT] 
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RESULTS AND DISCUSSION 



We have performed different calculations in order to 
extract the most relevant features in a cosmological N- 
body framework. We will first discuss the theoretical 
uncertainties on fluxes that we can infer from an analy- 
sis using WIMP fiducial models, combining the N-body 
data with analytical extrapolations. Afterwards, we will 
focus on the predictions that we can derive out of an av- 
eraged description of our 3D dark matter density map 
for different BSM dark matter particle candidates. 



A. Effects of the dark matter distribution and 
inhomogeneities 

1. Modifications of the spherical smooth profile 

Contrarily to dark matter contributions to the Galac- 
tic gamma-ray flux, it has long been shown that changes 
in the inner and outer logarithmic slopes of the smooth 
dark matter density profile have no significant effects 
on the charged cosmic ray signatures. This is due 
to diffusion, which tends to dilute the signal when 
integrated over large volumes (low/high energies for 
positrons/antiprotons). This has been checked for this 
N-body simulation, for which different fits of the den- 
sity profile give roughly the same statistical description 
of the data. This is illustrated in Fig. [3l where wc plot 
the positron and antiproton fluxes corresponding to three 
different fitted profiles, and where we also plot the result 
obtained with the common NFW model very often used 
in the literature to describe our Galaxy. The former three 
are characterized by an averaged density of pq = 0.25 
GcV/cm'^ at 8 kpc from the GC, while the latter is taken 
with its standard value of = 0.3 GeV/cm"^. It is note- 
worthy that there is absolutely no observable difference in 
the computed fluxes between the former three, and that 
the latter only gives a very small additional contribution 
at high (low) energy for positrons (antiprotons). This 
is merely due to the increase in the local density that 
translates to a factor of (0.3/0.25)^ to fluxes for small 
propagation lengths. This clearly demonstrates that even 
significant modifications in smooth profiles lead to minor 
effects, provided the local environment is not too much 
affected. 



2. Departure from spherical symmetry 

In the previous paragraph, wc have shown that the 
flux predictions were not sensitive to modifications of 
smooth spherical profiles. However, it turns out that 
though dark matter halos arc found to be roughly spher- 
ically symmetric in N-body simulations in average, they 
usually do not exhibit such a symmetry when 3-D scru- 
tinized. Furthermore, both observations and dyiiamical 
studies support axisymmctry instead [tO, [71|, . The 
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Figure 3: Comparison of positron and antiproton fluxes ob- 
tained with different smooth profiles: three are fitted on the 
HORIZON simulation (best fit, NFW and cored) and the 
fourth is the standard NFW fit very often used for the Milky 
Way. 



N-body framework allows to directly take the effects of 
departure from spherical symmetry into account. This 
may be seen in the left panel of Fig. 21 where we com- 
puted the positron and antiproton contributions associ- 
ated with the three projected density maps of Fig. [1] (an 
elliptic distortion appears in the xy plane, along the y 
axis). We see that for positrons (antiprotons respec- 
tively) the low (high) energy parts of the spectra are 
superimposed, whereas the high (low) energy parts dif- 
fer significantly. This comes from that low (high) en- 
ergy contributions correspond to a signal averaged over 
large volumes, erasing information on the source geome- 
try. Conversely, local variations of the dark matter den- 
sity are only relevant at high (low) energies for positrons 
(antiprotons). Therefore, global geometry modifications 
such as departure from sphericity can have an impor- 
tant impact if the averaged local dark matter density is 
changed accordingly. This can be further checked by ana- 
lytically deforming a spherical halo to model an elliptical 
one. The recipe is presented in [73| . where an elliptical 
profile of eccentricity a along the x axis in the xy plane 
can be inferred from the spherical dark matter density 
according to Pcmpticai(a;, y, 2;) = /Osphoricai(a;/a, ay, z). Of 
course, elliptically modifying the halo shape modifies the 
local dark matter density in the same time. Fixing the 
Earth position at x = 8 kpc, we have considered two 
cases: a = 2 and a = 1/2, which results in a modification 
of the local density by a factor of 3.26 and 0.26 respec- 
tively. Such modifications are expected to magnify the 



Benchmark point 


(GeV) 




{av){cm^^s ^) 


Branching ratios 


SUSY light- 1 : hght neutralino, annihilation through A funnel 
(scenario A \53} like) Non universal parameters 


9.7 


0.116 


2.17 X 10"^^ 


BRixX bb) = 0.873 
BR{xx ^ T+T-) = 0A25 


SUSY light-2 : light neutralino, annihilation into r"'"r~ via f 
(scenario B [53] like) Non universal parameters 


34.1 


0.129 


1.48 X 10"^*^ 


BR{xx^r+r~) = l 


SUSY-SPS2 : focus point, mixed higgsino-gaugino neutralino 
(SPS2-like) 

mo =2850 GeV, mj/a =500 GeV, Ao^O GeV, tan/3= 10 /i > 


198.7 


0.102 


2.30 X 10"^^ 


BRixx -> tt) = 0.719 
BRixX ZZ) = 0.0763 
BRixx -» ss) = 0.00387 


SUSY-SPS3 : neutralino-stau coannihilation 
(SPSS-like) 

mo =150 GeV, mi/a =680 GeV, Aq=0 GeV,tan/?= 10 /x > 


283.8 


0.102 


2.08 X 10"^^ 


BRixx -> bb) = 0.733 
BRixx ^r+T-)= 0.176 
BRixx -» = 0.0762 


SUSY-SPS4 : bino neutralino with main annihilation through Higgs funnel 
(SPS4-like) 

mo =770 GeV, mi/a =600 GeV, ^o=0 GeV, tan/3= 50 a* > 


252.2 


0.107 


1.50 X 10"^*^ 


BRixx bb) = 0.849 
BRixx ^ T+T-) =0.147 
BRixx -» ss) = 0.00387 


KK-1 : KK right-handed neutrino 

(warped-GUT-like). 

KK boson mass scale Mkk = 6 TeV 


50 


~ 0.1 


2.04 X 10"^^ 


BRixx bb) = 0.7 
BRixX l^l') = 0.035 X 3 
BRixx vv) = 0.065 X 3 


KK-2 : first excitation of the B*-^-* gauge boson 
(UED-like) 


1000 


~ 0.1 


1.7 X lO"^** 


Bi?(xX^ = 0.19x3 
BRixx -'{qq)u)= 0.1x2, 
BRixx vv) = 0.04 X 3 
BRixx ^ 0.01 X 3 


IDM-1 : light dark matter, annihilation through A funnel 

fi2 =35 GeV, AMao =10 GeV,AM^+=50 GeV, mh= 120 GeV, Aa = 0.1 


45 


0.109 


1.23 X 10"^'' 


BRixx -» bb) = 0.871 
Si?(xX^^^^") = 0.0435 
BRiHH cc) = 0.0835 


IDM-2 : heavy dark matter, annihilation into gauge bosons 

fi2 =1005 GeV, AMao =5 GeV,AM^+=10 GeV, mft= 120 GeV, Aa = 0.1 


1000 


0.117 


4.04 X 10"^^ 


BRiHH + W-) = 0.589 
BRiHH ZZ) = 0.223 
BRiHH hh) = 0.174 


LHM-1 : L- branch 

/ = 702 GeV, rriAH = 103 GeV, tuwh = 450 GeV, 
mz„ = 450 GeV, m* = 484 GeV, m^ = 120 GeV 


103 


~ 0.13 


1.4478 X 10"^'^ 


BRiHH -> W+W) = 0.7343 
BRixx -> ZZ) = 0.2657 


LHM-2 : U-branch 

/ = 1050 GeV, niAH = 162 GeV, = 678 GeV, 
mzH = 678 GeV, m* = 2410 GeV, = 400 GeV 


162 


~ 0.13 


1.9 X 10"^'' 


BRixx W+W-) = 0.6893 
BRixx -> ZZ) = 0.3107 



Table II; Benchmark points used in the calculation for Fig[7l 
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differences in fluxes observed in tfie left panel of Fig. |4l 
The full calculation indeed leads to the right panel of 
Fig. m where the former (latter) case is clearly shown to 
enhance (decrease respectively) the positron and antipro- 
ton fluxes by 1 order of magnitude, but mainly at energies 
corresponding to short propagation lengths. The effect 
is much less important at energies corresponding to large 
propagation lengths, where the spherical halo predictions 
are asymptotically recovered. The two previous examples 
are extreme configurations compared to the current ob- 
servational constraints, and are likely to be attenuated by 
the measured angle of ~ 20° between the GC - Earth 
axis and the semi-major axis [7(j|. 



3. Density fluctuations and sub-halos 

Beside global geometry considerations, dark matter in- 
homogeneities are expected to play a much more impor- 
tant role in characterizing the dark matter contribution 
to the antimatter spectra, but mostly at energies corre- 
sponding to small propagation scales. Indeed, at energies 
for which the propagation scale is large, the argument is 
the same as previously, and diffusion erases the effect 
of the source term fluctuations (contributions from the 
Galactic center dominate). On the contrary, the asymp- 
totic limit of fluxes for energies E corresponding to very 
short propagation scales only depends on the local dark 
matter density. As already stressed before, the main in- 
terest of using an N-body framework is to access a re- 
alistic inhomogeneous distribution of dark matter. By 
setting the terrestrial observer at different places in the 
N-body galaxy, we can survey the effect of density spatial 
fluctuations, that is the effect of modifying the local envi- 
ronment. This is illustrated in Fig. [51 where we plot the 
flux predictions for positrons (left panel) and antiprotons 
(right panel) associated with 100 different positions of the 
Earth within our N-body grid (the galactocentric radius 
is fixed at 8 kpc). For comparison, we report on the 
same figures the results when using the best-fit smooth 
profile, as well as the existing data, (from [tI . [tsI . [76j 
for positrons, and from [tI, [TJ, [zl, [zl, HO] for antipro- 
tons). What is interesting is, again, that contributions 
at energies corresponding to large propagation scales 
(high/low for antiprotons/positrons) are very well fea- 
tured by a smooth description, while fluxes corresponding 
to small propagation scales fluctuate a lot. Such fluctu- 
ations are due to modifications of the local dark matter 
density. The larger fluctuations come from configura- 
tions for which the local density is not constrained to be 
within [0.2, 0.4] GeV.cm^^^ (green curves - the black ones 
are density-constrained). It is easy to quantify the am- 
plitudes of these fluctuations, since the asymptotic flux 
(small propagation scales) merely scales like p©^. The 
value of pq = 0.25 GeV.cm"'^ used for the smooth profile 
allows to have an idea of the local dark matter densities 
corresponding to the extremal fluxes showed in Fig. [SI 
which spread over two orders of magnitude (a factor of 



2 for maximal, and a factor of 1/4 for minimal local 
densities). 

Although we can nicely deal with density fluctuations 
in our N-body galaxy, the lack of resolution at small 
scales prevent us to resolve sub-halos less massive than 
~ 10^ Mq. Moreover, more massive sub-halos are only 
found at the periphery of our simulation data, well be- 
yond the cosmic ray diffusion zone, and we do not have 
evidence for other sub-halos in the solar neighborhood 
that could affect the antimatter flux predictions. Sub- 
halos are expected to have much more important impacts 
on predictions because they are not merely density fluc- 
tuations: they are structured objects with inner profiles, 
the first virialized objects of the universe, of which the 
central density could be as or even more cuspy than the 
galactic center itself. The source term of antimatter pro- 
duction scaling like the squared density, sub-halos are 
therefore likely to constitute a major component of the 
exotic primary fluxes. Of course, some better simulations 
have already achieved resolutions good enough to resolve 
clumps as light as ~ lO^M© (e.g. the Via Lactea simula- 
tion described in [l^l), but they are still very far from the 
typical minimal masses ~ 10~®Mq set by the free stream- 
ing scales inferred from particle physics. Consequently, 
the intrinsic resolution limits of simulations make the use 
of analytical extrapolations unavoidable in order to study 
the effects of the smallest dark matter structures on pre- 
dictions. Besides and fortunately, the case of antimatter 
cosmic rays is well suited for such approaches, because 
contrarily to gamma-rays, the precise location and inter- 
nal features of any source are by no means observable 
with such messengers. Furthermore, an analytical and 
statistical recipe for calculating any clump configuration 
contribution to the antimatter primaries has been very 
well defined and described in [l^ and The latter 

reference provides an exhaustive and detailed method to 
deal with sub-halos as characterized in A— CDM cosmol- 
ogy. Armed with these analytical device, we can extrap- 
olate our results to different sub-halo configurations, all 
typified by a M~°''^' mass distribution, inner profiles and 
concentrations, and eventually by a spatial distribution. 

We have considered two specific cases. The first one 
is an optimistic maximal case in which clumps arc as 
light as 1O~^M0, with au = 2, inner NEW profiles, and 
obey the so-called Bullock - BOl - concentration-mass 
relation (8l| : in this configuration, we make sub-halos 
spatially track the smooth profile. The second one is 
approximately reminiscent from what the authors of [24 1 
found in their Via Lactea simulation, and will be referred 
to as the Via-Lactea-like configuration. It is featured by 
a minimal clump mass of IO^Mq, aM = 1.9 - i.e. ~ 10^^ 
times less clumps than in the maximal case inner NEW 
profiles, BOl concentrations, and a spatial distribution 
antibiased with respect to the smooth profile (see [s^ for 
more details on the spatial antibias). The reader is re- 
ferred to [20| for more insights on the impact of the clump 
parameters on predictions for the antimatter cosmic ray 
fluxes. We plot the positron and antiproton fluxes for 
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Figure 4: Positron and antiproton fluxes. Left panel: those corresponding to the three different Earth positions illustrated in 
Fig. [T] Right panel: taking into account an elliptical deformation of the dark matter halo (see the details in Sect. IV A 21 All 
fluxes are computed with the fiducial WIMP models (see Sect. II VA]) . 
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Figure 5: Varying the Earth position in the N-body simulation. 100 random positions are sorted and fluxes are computed for 
positrons (left) and antiprotons (right). Fiducial WIMP models detailed in Sect. IIV"^ are used. 



both configurations in Fig. [6l which somehow quantifies 
the consequences of our lack of spatial resolution, and 
might even aUow to by-pass it. These fluxes are char- 
acterized by mean values and associated statistical vari- 
ances (shown as contours) taking into account probabil- 
ities for clumps to be located inside the relevant (energy 



dependent) diffusion volume and to contribute to the sig- 
nals; in both configurations, the smooth halo is described 
by our best fit parameters (see Sect. |TT|. The maximal 
configuration (left panel) unveils a flux enhancement (so- 
called boost factor) with respect to the contribution of the 
smooth dark matter component which depends on the 
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energy, and saturates at energies corresponding to small 
propagation scales. Indeed, the smooth central halo con- 
tribution strongly dominates the signals against clumps 
at large propagation scales, whereas clumps are likely to 
overcome at small propagation scales if numerous enough 
in the local environment |l9l . [20l | . The maximal boost fac- 
tor is 4, which gives an idea of what we could expect 
from such a small scale dumpiness. Taking a r~^'^ pro- 
file inside clumps instead of NFW would provide an addi- 
tional factor of ~ 10. The Via-Lactea-like configuration 
does not provide any flux enhancement in average, but 
the statistical variance of signals is still relevant at short 
propagation scales, for which a very close object, a very 
improbable situation, could increase the signals. Notice 
that the sub-halo effects yield uncertainties of the same 
order of magnitude as the density fluctuations and de- 
parture from spherical symmetry previously considered. 



B. Predictions for different dark matter particle 
candidates 

After having tried to quantify the effects coming 
from possible dark matter density fluctuations with a 
WIMP model-independent setting, it is useful to restore 
those predictions in a more constrained particle physics 
scheme. This really allows to define and compare the 
flux orders of magnitude that WIMP searches do have to 
deal with. To this aim, we have used an Earth location 
in the N-body frame which provides intermediate fluxes 
as illustrated in Fig. [5l with a local density of pq ~ 0.25 
GeV.cm-3. 

We have derived predictions for all dark matter can- 
didates presented in Sect. IIVI of which the features are 
summarized in Tab. |TT] - some of them are detectable 
at the LHC, with well defined related signatures. The 
positron (left) and antiproton (right) fluxes are reported 
in Fig. El for SUSY (top) and non-SUSY particles (bot- 
tom panels, respectively). On the sa me p anels are pic- 
tured the existing data (from [tI . [tsI . [76| for positrons, 
and from [tgI. FttI. TtsI. [tqI. ISOj for antiprotons). In the top 
panels, we also report the secondary backgrounds com- 
puted by [131 for positrons (top left panel), and by [11] 
for antiprotons (top right panel), both using the same 
median set of propagation parameters employed here for 
primaries. From this figure, we clearly see that only the 
lightest SUSY candidate is about to imprint the antipro- 
ton spectrum, but at very low energies (< 1 GeV) , a spec- 
tral region which very much suffers from solar modulation 
effects (the scatter in the low energy BESS data points is 
mainly due to different solar activities during the obser- 
vations). All other candidates remain far lower than the 
existing data, and are therefore poorly constrained. Even 
a one or two order of magnitude enhancement, due to the 
various effects considered above, would not be enough 
to detect nor constrain the studied dark matter candi- 
dates with the current data. We finally stress that even 
primary signals high enough to be detectable would be 



constrained to not exceed too much the secondary con- 
tributions, which are already in rather good agreement 
with the data. 

Another trendy quantity which has recently strongly 
drawn attention is the positron fraction, that is the ratio 
0e+/(0e- + 4'e+) — whcrc (/)g- is the electron flux. The 
focus has hold since August (2008), when the PAMELA 
collaboration has presented some preliminary measure- 
ments up to 270 GeV at the IDM08 conference in Stock- 
holm, which indicate that the positron fraction signif- 
icantly increases with energy contrarily to what is ex- 
pected e.g. from the predictions of [83| of both the as- 
trophysical electron and secondary positron cosmic ray 
fluxes. These preliminary results confirm the trend pre- 
viously observed e.g. in the HEAT and AMS experi- 
ments [13, mi, nil. Furthermore, newer predictions of the 
secondary positron flux have just been published in [13, 
which have been derived in the same slab diffusion model 
as we use in this paper. The authors have also calculated 
the theoretical uncertainties coming from various origins, 
and found a rather good agreement between their pre- 
dictions of the secondary positron flux and the existing 
data on the positron flux (see Fig. [71 top left panel). In 
their discussion on the positron fraction, they have em- 
phasized the lack of knowledge of the electron cosmic ray 
flux above a few tens of GeV, where it is not well con- 
strained. They have nicely shown how. when describing 
the high energy part of the electron flux with a single 
power law, very slight modifications of the spectral index 
can lead to quite different interpretations of the positron 
fraction data, spanning from a large excess to no excess 
at all with respect to the predictions. Though we have al- 
ready shown in Fig.[7]that our predictions for the primary 
positrons lie well below the data and the expected back- 
ground, we also plot in Fig. [5] the corresponding positron 
fraction for the sake of completeness. To this aim. we 
have modeled the electron flux like in [13, with a fit on 
the AMS data and constrain the spectral part above a 
few GeV to be a single power law with index -3.44 corre- 
sponding to the best fit value obtained in (STj . in which 
data up to ~ 1 TeV were utilized. Note that this electron 
spectrum differs from that of [27j , since it does not corre- 
spond to their hard, neither soft indices, but to a median 
one instead. In Fig. [51 w e also report the positron frac- 
tion data released in [83, [H, [s^ , and we also show the 
PAMELA data but up to 10 GeV only [ll], since higher 
energy points are not yet officially released. The positron 
fractions associated with the SUSY dark matter models 
are plotted in the left panel of Fig. [51 those correspond- 
ing to the non-SUSY models are featured in the right 
panel. We see again that all primary signals are orders 
of magnitude far below not only from the data but also 
from the background expectation. Modifying the elec- 
tron spectrum would not change our conclusions here. If 
one still tries to accomodate the positron fraction fiux by 
tuning the dark matter properties, the counter part in 
antiprotons will be a strong limit, unless the antiproton 
production is suppressed, which is the case for our light 
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Figure 6; Analytical predictions for the positron and antiproton fluxes when considering sub-halos in the Horizon framework. 
Left panel: maximal configuration, where 10~^Mq NFW clumps are involved. Right panel: Via-Lactea-like configuration, 
where we consider IO^A/q NFW sub-halos. Fiducial WIMP models detailed in Sect. llVAl are used. 



neutralino model of 35 GeV for which only the Icptonic 
couplings are relevant (see Tab. |TT|. Such a conclusion 
has also been stressed in (soj . 

All statements in this subsection rely of course on the 
fact that the propagation model we used is defined by 
the median set of propagation parameters of Tab. [J A 
more favorable set, which for instance would involve a 
larger diffusion coefficient K and a larger vertical exten- 
sion L of the diffusion zone (such that K/L, the main 
constrained quantity, remains constant), could lead up to 
an additional order of magnitude in terms of fluxes (sec 
e.g. [13, HH, nil ) I ^^'^ could therefore make of the current 
data very strong limits. This points towards the neces- 
sity to have more precise and over larger energy ranges 
data of secondary/primary cosmic rays to much better 
bound the propagation allowed configurations. Further- 
more, we see on Fig.[7]that there is an overwhelming lack 
of high energy data above ~ 10 GeV, precisely in spec- 
tral regions where the astrophysical background coming 
from secondaries is expected to be power law suppressed. 
Higher energy data could still exhibit spectral transitions 
due to some exotic components, provided the dark mat- 
ter candidate is massive enough. However, this would re- 
quire dark matter particles and/or spatial distributions 
with quite different properties from those used in this pa- 
per, since we have have shown them tho have contribu- 
tion up to several orders of magnitude below the existing 
data. The PAMELA satellite is currently improving the 
measurements on antiprotons and positrons in a slightly 
higher energy range than before [901, and the GLAST 
satellite will also have the potential to detect high en- 



ergy electrons and positrons [9l| 



C. Limits of the N-body framework 

N-body simulations are the only way to scrutinize 
structures that have virialized in the non linear regime, 
and especially when one is interested to check some pe- 
culiar effects, like those coming from inhomogeneities or 
asphericity. Nevertheless, such a framework still suffers 
from many drawbacks, due to their incompleteness. First 
of all, the limited spatial resolution prevents the study 
of the smallest scales that substructures can have, and 
which are fixed by the free streaming length of dark mat- 
ter particles in the early universe. Some authors have 
tried to tackle this issue, by performing simulations re- 
solving ^ IO^^Mq objects, but their results hold only 
at a redshift z — 26, when galactic halos are not yet 
formed [o^- Even by taking the fastest available su- 
percomputers for cosmological simulations, the authors 
of [il] only reached a resolution allowing the identifica- 
tion of ^ IO^Mq objects, that is more than 10 orders 
of magnitude beyond what is expected from WIMP cold 
dark matter. 

Moreover, it is noteworthy to recall that baryons are 
likely to play an important role in the dynamics of galax- 
ies, and they are not yet fully incorporated in cosmolog- 
ical simulations down to the galactic scale. For instance, 
galactic bars observed in the centers of spirals have long 
been shown to exchange momentum with the dark mat- 
ter halo, and are therefore expected to modify the dark 
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Figure 7: Left: positron fluxes for SUSY dark matter (top) and non-SUSY dark matter (bottom). Riglit: same as left panels 
but for antiprotons. For comparison, the secondary background for positrons (antiprotons, respectively) as calculated by [23] 
([2^) is reported on the top left (right) panel. 



matter distribution (see e.g. a recent analysis in [73 |. 
and references therein). Moreover, it is well known, but 
scarcely mentioned in the present context, that the avail- 
able data on the dynamics of the Milky Way, the target of 
interest in this paper, strongly disfavor cuspy halo models 
(see more detailed arguments in e.g. [93l|). 

In order to illustrate this latter statement, we have 
used the kinematic data provided by [qJ], for which 
the authors subtracted the contribution of the baryons 
as described by the bulge-disc-bar model of [7l|. In 
Fig. [S] (left panel), we have reported the circular ve- 



locities against the galactoccntric radius, inferred from 
the smooth dark matter profiles fitted on our N-body 
data and presented in Sect. [TTl from the standard NFW 
profile, and also from a ad hoc profile with parameters 
a = 2, ^ = 1, 7 = 0.5, r, = lOkpc, pq = O.SGeV.cm-^ 
(cf. Eq.m. The latter has no physical justification at all, 
and has only been tuned to correctly adjust the velocity 
data. For comparison, we have plotted the same pro- 
files against the density data extracted from our N-body 
simulation (right panel). 

Disregarding the ad hoc profile for the moment, it is in- 
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Figure 8: Positron fraction for SUSY dark matter (left panel) and non-SUSY dark matter (right panel). The secondary positron 
background is that calculated in [23| with the median set of propagation parameters. 



teresting to remark that while the standard NFW hardly 
better fits the velocity curves than the others at large 
radii {r > Q kpc), all the cuspy profiles fail to repro- 
duce the dark matter rotation curve in the inner 4 kpc 
— the central parts of our Galaxy are influenced by the 
bar, which perturbs the kinematics of the gas; thus the 
points inside about 3 kpc are not to be considered. These 
cusps systematically lead to an overestimate of the dark 
matter mass. Between 4 and 8 kpc the fits are also bad, 
but take advantage of the large dispersion of the observa- 
tional points to remain within the cloud of points. How- 
ever, they adjust the highest values at the inner radii 
and the lowest values at the largest radii, i.e. are not 
acceptable in the sense. Beyond 8 kpc the values pre- 
dicted by the cuspy profiles are by far too low compared 
to the observations, even when one takes into account the 
considerable uncertainties of the observations outside the 
solar radius. The ad hoc profile, which is close to a cored 
isothermal profile, provides a good fit to the observations 
though it may overestimate the velocities at large radius 
(see Fig.O left). This is not surprising since it was tuned 
for this purpose. Nevertheless, interestingly, this profile 
does not fit the N-body data (see Fig. [9l right). This 
underlines one of the well known and much discussed 
discrepancies between observations and the cosmological 
N-body simulations in the scale of galaxies (see e.g. [IHl 
and references therein). Such discrepancies are believed 
to be due to the fact that cosmological N-body simu- 
lations do not include yet the baryonic component, of 
which the presence could lower the cusp into a core. 

Hence, we stress that the predictions we provided in 
the previous section should be better considered as in- 



dicative, because they do certainly not contain the full 
theoretical uncertainties, especially those coming from 
our ignorance of the baryon influence. Nevertheless, our 
calculations, as performed in an N-body context, still al- 
low a better understanding of the uncertainties which 
come from the possible fluctuations of our local dark 
environment. Computing the fluxes by using the pre- 
viously discussed ad hoc proflle would give exactly the 
same trend as those reported in Fig. [31 with a slightly 
increased asymptotic value at short cosmic ray propaga- 
tion scales, by a factor of {pf^°'' / p]^°"^°''f = 4. In or- 
der to not further complexify the reading, we do not plot 
the corresponding predictions, since this profile does not 
rely on any theoretical motivation. Eventually, note that 
despite the already mentionned drawbacks, our analysis 
sketches the bases for further studies using more compli- 
cated frameworks (including baryons for instance). 



VI. CONCLUSION 

In this paper, we have studied the theoretical un- 
certainties associated with dark matter density fiuctu- 
ations affecting the primary antimatter cosmic ray fluxes 
(positrons and antiprotons), possibly produced by dark 
matter annihilation in the Galaxy, by using an N-body 
simulation together with analytical descriptions of a 
galaxy-like dark matter halo. This is the first attempt 
in this research field to directly connect the source term 
of the cosmic ray propagation equation to a 3D den- 
sity map coming from a cosmological N-body framework, 
while this has already often been used for predictions 
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Figure 9: Left; Radial velocities inferred from the spherical dark matter halo models discussed in Sect.]!!] plus an ad hoc profile 
with a = 2, /3 — 1, J — 0.5, Ts = lOkpc, pq = 0.5 GeV.cm"'^. The black stars represent the kinematic data from which the 
baryonic contribution has been subtracted; cf. Fig. 3 of [oj. Right: same density profiles reported on the N-body density 
data, for comparison. 



of gamma-ray fluxes. Indeed, the latter situation does 
not involve cosmic ray propagation, which seriously com- 
plicates the calculations, and makes them much more 
time consuming. The main purposes of this work was to 
quantify those uncertainties, and in some cases to try to 
understand them by analytical means. While this was 
achieved in a WIMP model-independent setting, we also 
aimed at reviewing the status of predictions for many dif- 
ferent well motivated dark matter candidates in particle 
physics, supcrsymmetric or not. 

The N-body framework led us to study several ef- 
fects, such as non-spherical dark matter distributions 
(Sect. IVA2|) and density fluctuations (Sect. IV A3p . 
which are not accounted for when using the traditional 
spherical smooth halo models. We have shown that, as 
expected, fluxes are mainly affected at cosmic ray en- 
ergies corresponding to short propagation characteris- 
tic lengths (low/high for antiprotons/positrons) and to 
local contributions, while when looking at energies of 
large propagation lengths, the smoothing due to diffusion 
erases the peculiarities in the spectra, and is well repro- 
duced by a spherical halo model. The main uncertainties 
on the predicted fluxes come therefore from those on the 
local dark matter environment, which we have shown to 
fluctuate a lot in the N-body framework. This can lead 
to ± 1 order of magnitude in terms of flux, at high energy 
for positrons, and low energy for antiprotons. This could 
be more constrained with much better limits on the local 
dark matter density. 

Moreover, though our simulation has a too poor spa- 



tial resolution to resolve substructures inside the galactic 
halo, we have extrapolated our results in order to include 
the potential effects of the presence of sub-halos by using 
the analytical method detailed in [2^. We have consid- 
ered two situations: one maximal, analytically extending 
the clump mass spectrum down to 10~^ Mq, led to a 
flux enhancement by an energy-dependent factor reach- 
ing 4 at maximum; a second, Via-Lactea-like, involved 
clumps down to IO^A/q only and did not result in any 
flux enhancement. 

Consequently to the previous points, the smooth and 
spherical modeling of the dark matter halo leads to a very 
good approximation of a more complete and complex sit- 
uation, as derived from an N-body simulation, at large 
characteristic propagation lengths for antimatter cosmic 
rays (low/high energies for positrons/antiprotons). Nev- 
ertheless, though predictions in the spectral regions as- 
sociated with those large characteristic scales are also 
the less affected by the uncertainties coming from dark 
matter density fluctuations, we remind that they are still 
very sensitive to the used propagation model. 

Beside this trial for directly measuring the theoret- 
ical uncertainties with an N-body experiment, we re- 
viewed the predictions of the positron and antiproton pri- 
mary contribution associated with some specific particle 
physics dark matter candidates (Sect. rVBj) . We have not 
only included some popular supersymmetric and extra- 
dimensional models, but also some more specific ones 
based on minimality arguments, such as the little Higgs 
model, or the inert doublet model. We have shown that 
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predictions are well lower than the existing data, except 
perhaps for the antiproton flux associated with the light- 
est neutralino model, just because of the more favorable 
1/m^^ factor appearing in the flux expression. Neverthe- 
less, even if the latter appears in tension with the cur- 
rent data, its contribution to the antiproton flux occurs 
at energies where solar modulation effects are expected 
to be important, and it would be hard to find a clear 
signature in this spectral region. The other models are 
well below the experimental measures, but higher energy 
measurements could perhaps unveil some spectral transi- 
tions followed by an energy cut-off which could possibly 
be attributed to dark matter. In this case, uncertain- 
ties associated with the dark matter density fluctuations 
would be much less stringent for the antiproton signal 
than for the positron because of propagation scale argu- 
ments. On the contrary, a nearby clump would favor a 
detection with positrons at high energy, much more than 
with antiprotons. Anyway, higher energy data are neces- 
sary in this field, for these antimatter species as well as for 
standard nuclei species which will yield much better con- 
straints to the propagation modeling. Obviously, indica- 
tions of a WIMP candidate mass, as expected from the 
LHC, would provide the relevant spectral region where to 
concentrate the astrophysical searches, and would allow 
to concentrate on much more subtle effects. 

Finally, we made a self-criticism exercise - Sect. IV Cl - 



by comparing our different halo models to the star kine- 
matic data available for our Milky- Way Galaxy, after sub- 
traction of the baryonic contribution as modeled by [o^l ■ 
We show that the dark matter contribution to the veloc- 
ity field as inferred from our N-body simulation, as well 
as more general ones, are far to reproduce the kinematic 
data, systematically leading to a mass excess in the cen- 
tral regions (< 4 kpc) of the Galaxy. This clearly calls 
for restraint in our potential claims, and clearly points 
towards some large amount of ignorance in the intimate 
(and gravitational) relations between dark matter and 
baryons. 
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